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Abstract 

Using a pair of the lightest mirror nuclei, 3 He and 3 H, we study the effect 
of the medium modification of pion fields on the flavor non-singlet structure 
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The partonic distribution functions of the nucleon, in particular the flavor dependence 
of the antiquark distributions, are of considerable interest [Q]. Within the framework of 
perturbative QCD (pQCD) the light quark sea is expected to be flavor symmetric. However, 
the experimental data || contradict this idea, revealing an excess of d over u in the free 
proton. This inconsistency indicates that non-perturbative effects should be responsible 
for the flavor asymmetry in the light sea quark distributions. For example, some flavor 
asymmetry was anticipated before the measurements on the basis of the chiral structure of 
the nucleon 0. 

The physical proton has a relatively large 7r + -neutron Fock component which naturally 
leads to a surplus of d It is known that this Fock component offers the main contri- 

bution to the d excess and that the contributions of the other mesons and A isobars have 
opposite signs and tend to cancel each other Jl[. An alternative explanation for an excess 
of d over u involved the Pauli exclusion principle, given that there are two valence u quarks 
in the proton and one valence d [|J. Perturbative estimates failed to support this ||, with 
the first non-perturbative explanation of the origin of such an effect in terms of the vacuum 
structure of the proton given by Signal and Thomas ||. Estimates of this effect within chiral 
quark models have also been given in Refs. f?]-|9| • ^ may well be that the experimentally 
observed excess involves contributions from both of these effects [|TIJ . 

One way to learn more about the non-perturbative structure of the nucleon is to study 
the non-singlet difference between the proton (p) and neutron (n) structure functions, for 
nucleons bound in a pair of mirror nuclei [jTT] , |T2|| - In this case any discrepancy between 
theoretical predictions and observed data will indicate a modification of the non-perturbative 
mechanism giving rise to the flavor asymmetry in the free proton, in the nuclear medium. 
In particular, such a discrepancy would be a sensitive probe to study pions in nuclei. In this 
paper we examine the effect of changes in the pion cloud on the non-singlet combination of 
nuclear structure functions, using the lightest pair of mirror nuclei, 3 He and 3 H. 

How is the pion field modified in a nucleus ? To study it we concentrate here on only Fock 
states consisting of a "bare" nucleon and pion, and ignore nuclear binding, Fermi motion and 
shadowing/ ant ishadowing effects for the moment. Under these assumptions, the structure 
functions of the proton and neutron in the nucleus A are given by [I[ 

F<l = Z p/AF<2 + fnOp/p/A ® F% + fn°p/p/A ® F£ + f-^+n/p/A ® +n/p/A 

®Ff, (1) 

F^ A = Z h /aF2 + f n o n / n /A <8> F£ + f^On/n/A <g> F% + f n - p / n /A <£> F% + fn-p/n/A <£> F% , (2) 

where F% is the structure function of a 'bare' proton (neutron). The probability to find the 
'bare' proton (neutron) in the physical proton (neutron) in A is denoted by the normalization 
constant, z p ( n )/A- The shorthand notation, fuB/N/A <8> F 2 K , stands for the convolution of the 
(light-cone) momentum distribution of the pion M (= 7r~, 7r°, 7r + ) per N (= p, n), fMB/N/A{y) 
(B = p,n), and the structure function of K, F^{x) (K = B, M) [0: 

rl—% f X \ 

f mb in I a ® Ff(x) = ^ dyf M B/N/A(y)F2 f -^—^ J , (3) 

fus/N/A ® F 2 M (x) = £ dyf MB/N/A (y)F 2 M (J j . (4) 

The nuclear structure function is then simply given by F 2 A (x) = ZF% (x) + NF£' A (x), 
where Z and iV are the numbers of protons and neutrons, respectively. 
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We consider a pair of mirror nuclei: A = Z + N [Z > N) (proton rich) and A' = 
Z' + N' (N f > Z') (neutron rich). In a nucleus we can expect a significant difference in 
the positive and negative pion light-cone momentum distributions. Recently, Korpa and 
Dieperink have calculated the pion fields in asymmetric nuclear matter ||13|| , with a result 



which is consistent with the Drell-Yan experiment of Aide et al. ||14|| . Their result suggests 
that the difference in the distributions basically comes from two factors. One is the Pauli 
blocking of the nucleon in the final state; in the proton rich nucleus A which we consider 
here, the emission of tt~ (from a neutron, creating a proton in the final state) is more 
suppressed than tt + emission. The other effect is the dressing of the pion propagator in 



matter ||15|| , where the particle-hole self-energy dominates. Korpa and Dieperink find that 
the delta-hole contribution is minor and that the neutral pion field is not much altered in 
the nuclear medium [[13]. In summary, their analysis suggests that in the proton rich nucleus 
A the 7r + (7r - ) field is enhanced (reduced) as compared with that in the free nucleon, while 
the 7T° field is not changed a great deal. 

In a nucleus the Coulomb interaction may affect the shape of the pion momentum dis- 
tribution and, of course, one cannot use isospin symmetry: 

fn+n/p/A 7^ fn~p/n/A' and f n -p/n/A 7^ fir+ n/p/A 1 ■ (5) 

However, for the reasons discussed above, we suppose that the 7r° distribution is not changed 
much in matter: 

f-n-Op/p/A = fn°n/n/A = fn°p/p/A' = fn°n/n/A' = fn°N- (6) 

Here f n o N is the tt° distribution in the free nucleon, which is given by the Sullivan process 



WvHt^tt; \o / dti ,^ \M + y*M$ r ), (7) 



9 2 r dh 2 gjv(g) 

167r 2 y(l - yf Jo '(M^-s) 2 ' 

with g(= 13) the it — N coupling constant, kf the transverse momentum squared of the pion 
and 

s = ml + kj + Ml + k\ 

y 1-2/ 

The free nucleon mass is denoted Mn (0.94 GeV) and m n (0.138 GeV) is the pion mass. 
The form factor, F n N(s), is given by |L| 



Fttn(s) = exp 



Ml 



2A 2 



(9) 



with A the cut off parameter. 

We divide the pion distribution into two pieces: 



f-K+n/p/A{A'){y) — fw+n/piy) + ^f^+/A{A') {v) j (io) 
U-p/n/A(A')(y) = U-p/n{y) + $fn-/A(A>){y), (ii) 

where f n +n/p(n-p/n) is the momentum distribution of 7i + (n~) in the free proton (neutron). 
The nuclear many-body effects on the pion field in A(A') are expressed by 5/m/a(A')- 



3 



The normalization constants in Eqs.(P and (0) can be related to those for the free 
nucleon. For example, using Eqs.(H) and fllOD we find 



Zp/A = 1 - (U°n) - (fw+n/p) - (SU+/A) = ZN - ($fn+/A), (12) 

where zn (= 1 — (f n °N) ~ (fn+n/p)) is the normalization constant for the free nucleon. Finally 
then the nucleon structure functions in those nuclei become: 

F V/A = p p _ {6U+/A) PP + §U+/A g, { pn + ^ (13) 

F 2 nM = F 2 " - (8U- /A )F? + 8 fa- /a ® [F p 2 + Ff], (14) 
Fl IA ' = Fl - {8U +/A ,)Fi + 8U +/A , ® [F? + Ff], (15) 
F 2 n/A ' = F 2 " - (5fr/A>)F? + 8fn-/A' ® [^ 2 P + ^1, (16) 

where F 2 p( ' Tl ' ) is the free proton (neutron) structure function. 

Now we study the lightest mirror nuclei: A = 3 He and A' = 3 H. In 3 He, the ir + meson is 
generated from the proton and the final state is given by 3 He = 2p + n — > p + 2n + 7T + , where 
7r + feels a repulsive force from the Coulomb interaction with the single p. (We neglect two 
pion emissions in the final states, such as 2p + n — > 3n + 2tt + , for which the probability is 
expected to be very small.) On the other hand, the ir~ meson is produced by the neutron, 
and the final state is 3p + it~, where the ir~ feels a strong attractive force due to the 3p — tt~ 
interaction. Thus, we expect that the Coulomb force between 3p and 7r~ is about three times 
stronger than that between p and ir + in the former case. In 3 H, the n~ feels an attractive 
force in the final state 2p + n + 7r~, and the Coulomb force is twice as large as that for the 
7r + in 3 He. On the other hand, the 7r + in 3 H does not feel any Coulomb force because the 
final state consists of 3n + 7r + . 

In order to evaluate Eqs. (p!3D -(|l6|), we need to estimate the distributions, 5/m/a(?/); m 



3 He and 3 H, individually. The Coulomb force may change the shape of the pion momentum 
distribution. As discussed above, since it acts on n~ as an attractive force in the nucleus, 
the wave function of the pion in coordinate space shrinks. This means that the pion gets a 
(relatively) higher momentum and the shape of the distribution in momentum space should 
shift toward larger y. For the 7r + the distribution should be modified the opposite way, 
because it feels a repulsive force. To calculate these effects quantitatively requires very 
complicated many-body calculations, including Coulomb forces. This is extremely difficult 
and in order to make a first estimate of the effects one might expect the following simple 
scaling assumption to be reasonable. That is, the change in the pion distributions are 
assumed to be given by 

8U+/A(y) = 2o w +/a(1 + /5 7 t+/a)/7t07v((1 + fl-K+/A)y) = 2aL ir +/Afi+/A(y)> (I 7 ) 

8fn~/A(y) = 2^-/^(1 - /VaO/7t°Jv((1 - Pk-/a)v) = 2 a Tx-/Afl-/A{y)> (i 8 ) 

8f n +/A'(y) = 2a 7r +/ J 4'/ 7r ov(2/), (19) 

8fn-/A'(y) = 2a w -/ A '(l - Ar-A4')/-7r°iv((l ~ P*-/A')v) = 2<**-/A'fZ-/A'(v)> ( 20 ) 

where aM/A(A') represents a change caused by the strong interaction in a nucleus (for ex- 
ample, Pauli blocking, correlations of random phase approximation (RPA), etc ||13| ) and 
{3m/a(a>)(> 0) describes a shift of the distribution because of the Coulomb force. Note 
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that the 7r + in 3 H does not feel any Coulomb force. As pointed out above, we expect that 
Ar+ 1 A '■ Pn-/A '■ Pn-/A' = 1:3:2 from the point of view of the strength of the Coulomb force 
acting on the pion. We therefore choose 3ft n +/A = ftn-/A — §Ar-/A' = ft > (ft is assumed 
to be small). The new function, fj^/A(A'Y * s normanze d as 



1 

d yfli/A(A')(y) = (/ttojv), (21) 

where we ignored a tiny quantity stemming from J x _g ; dyf n o N (y). 

Next, we suppose that the average number of pions per nucleon in a nucleus is equal 
to that in the free nucleon - experimental indications are that the pion field is not much 
enhanced in a nucleus The requirement of pion number conservation reduces the 



number of parameters "a" in Eqs. ([T7|)-(|20|) . We find that in 3 He, 2a w +/A + ol-^-ia = 0> 
while in 3 H, a^+/A> + 2a 7r -/ ( 4/ = 0. Thus, we set 2a w +/A = —ol^-ia = &a > (the n~ 
field is suppressed in 3 He) and 2a n ~/A> = —ol^+ia 1 = &A' > (the n + field is suppressed in 
3 H). Furthermore, since qca(A') describes the change of the pion field because of the strong 
interaction, we can set a a = &A 1 = oi (isospin is a good symmetry in this case). This leaves 
just two parameters, a and ft. 

We should note here that even if a = the Coulomb effect would modify the proton 
and neutron structure functions in the nucleus. Such a case could be described by replacing 
the pion distribution in the free nucleon structure function, /mb/n, in Eqs.(|I3|)-(|rfi|) with 
flt/AiA')- However, we expect that by itself the Coulomb effect on the structure function 
should be quite small (see below Eq.(|50D and Fig. |^). We therefore neglect the Coulomb 
effect on Ff W in F^ n)/A{A,) . Equations (jTH)-© then give 



5F V/ He _ F p/ He _ ^ = ^f^V + ^ ^ g, { pn + ^ 

7 s Be 



5Ff He = Ff He - = 2a(f n o N )F? - 2a£_ / , ire ® [Ff + Ff], (23) 



6f p/ h = F p/ h_ f p = 2a(f 7r o N )Fi - 2a/ 7r o JV ® [Ff + Ff], (24) 
5Ff H = Ff H - F% = -a(f n o N )F? + af* n - /3H ® [Ff + Ff), (25) 



where, 



f: +/ >He(y) = (i + \ft) U°n ((i + \ft) y) , (26) 

K-/>hM = (1 - /?) W(l - #v), (27) 

K-My) = (i - /.oiv ((i - ^) y) ■ (28) 

In Fig. p] the pion distributions provided by the scaling assumption are presented taking 
A = 1 GeV (see also below Eq. (j36j) ). As an example, we choose ft = 0.1, which means 
that, for instance, the wave function of the tc~ in 3 He shrinks by about 10% in coordinate 
space because of the Coulomb force. The negative pion distribution carries somewhat higher 
momentum, while the positive one shifts toward lower y, compared with f n o N . Taking the 
non-singlet combination of the structure functions of 3 He (F 2 He ) and 3 H (F 2 H ) we find 
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F^ - F 2 H = {Fl - if?) - Aa(f n o N )6F 2 N 

- 2a[/;_ /3H ® + /;- /3ffe ® *T - f: +/3He ® ^ + - ® ^1, (29) 

where 

= F 2 p (x) - i>(x) = l -x[u v {x) - d v (x)}, (30) 

with u v (d v ) the valence u(d) distribution in the bare proton. 

Figure |] illustrates the nuclear and Coulomb effects on the non-singlet structure function 
of the A = 3 system, which is given by 5F 2 A=3 = (F 2 He — F 2 H ) — {F 2 — F 2 ) . For the numerical 
calculations we have chosen (at Q 2 = 4 GeV 2 ) fTB] 



xu v (x) = 0.65452 x x a38 (l - a;) 2 ' 49 (l + 10.5a;), (31) 
xd v {x) = 0.028660 x x a07 (l - x) 4 ' 63 (l + 150x), (32) 

and F? + (x) = Ff{x) = 0.98863 x x a61 (l - x) im 0. Clearly the effect of the Coulomb 
distortion is quite small in the region x > 10~ 4 , even if we choose (3 = 0.2. (We have checked 
that the Coulomb effect on 5F 2 ^ A ^ A ^ individually is also small.) 
The nuclear Gottfried integral, Iq A (z), is defined by JTT 



I A < A \z) = U A -[F 2 A {x) - Ff\x)\ = I»{z) + 8I A ' A \z), (33) 

with Y{= Z — N) the number of excess protons in A and Iq(z) the Gottfried integral for 
the free nucleon. The nuclear effect is described by the second term on the r.h.s of Eq. (|33|) 

6l£ A '(z)= [ A ^8F A ' A '(x), (34) 

J z X 

where 5F A ' A ' = ^(F 2 A - F 2 A ') - (if - F 2 n ). In the case of the A = 3 system 5F A ' A ' is given 
by 5F 2 =3 . As we have already discussed in Ref. the Gottfried integral is generally 



divergent when the effect of charge symmetry breaking is included, even for the free proton 
and neutron [ |T2]j . 

If we set (3 = 0, 5F 2 =3 reads 

SF 2 A = 3 = -4a [(U N ) + U N ®} 5F 2 N . (35) 

Since the effects of the nuclear binding and shadowing are ignored in the calculation for the 
time being, the change in the Gottfried integral is convergent and it is given by 

sf G He:iH (0) = - 8 -a(Uo N ). (36) 

With the cut-off-mass A = 1 GeV in f w o N , we find (/^ojv) = 0.083. The Gottfried integral for 
the free nucleon is then estimated to be Iq(0) = |(1 — 4(/ 7r oj V )) = 0.223, which is consistent 
with the measured value (0.235±0.026) |l[. We can see that the modification of the pion field 
enhances the flavor asymmetry in a pair of mirror nuclei and hence it reduces the Gottfried 
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integral. If a = 0.05(0.1)[0.2] (and p = 0), 5f G He ' 3H {0) = -0.0111(-0.0221)[-0.0443], 
corresponding to a reduction of the Gottfried sum by about 5(9) [20]% from the free value. 

Next we consider the nuclear binding and shadowing effects. (For recent reviews see 
Ref. JIB] .) In Ref. [0 we studied shadowing and anti-shadowing corrections to the flavor 
non-singlet structure function using the Gribov-Glauber multiple scattering formalism. We 
found that the non-singlet structure function is enhanced at small x by nuclear shadowing, 
increasing the nuclear Gottfried integral (z is chosen to be 10 -4 in Eq.([3^)) by between 15 
and 41%. The enhancement of the non-singlet structure function is caused by the difference 
between the density distributions of 3 He and 3 H. 

In the shadowing region the structure functions of 3 He and 3 H are given by ]12 



F 3 2 He = 2F% + F 2 n - (2.5/ 3He - g 3He )Fg - 0.5fs He F^ (37) 
Fl H = Fi + 2F? - (2.5/sh - g* H )F? - 0.5f. H Fl (38) 

where fzHe( A H) and gsne( s H), respectively, describe the single and double rescattering pro- 
cesses in 3 He ( 3 H), which depend on the nuclear density distribution. It is also necessary to 
include anti-shadowing in order to produce the structure function around x ~ 0.1 (using the 



baryon number and momentum sum rules). Detailed discussions can be found in Ref. |fL2j| . 

Replacing the structure functions of the free proton and neutron in Eqs. fl3"7j) and (3S) 
by F^ A( ^ A > (given by Eqs.(]T3|)-(]T6D), we can calculate the non-singlet structure function of 
3 He and 3 H, including both the modification of the pion fields and nuclear shadowing. (This 
means that the effect of the pion cloud modifies the nucleon sea quarks only.) Note that this 
simple replacement is an approximation made in order to see the effect of the change of the 
pion fields. To treat this problem rigorously it would be necessary to construct a model where 



pions are handled consistently and contribute to both shadowing and anti-shadowing |L9 
which goes beyond the scope of the present work. 

In the calculation of nuclear shadowing, the ground-state wave functions of 3 He and 3 H 



are assumed to be given by gaussian functions in coordinate space ]l2j, |\l/| 2 oc exp[— r 2 / (26)], 
where the parameter b determines the correct matter radius of the nucleus. We take b = 
40.59(30.06) GeV" 2 , which produces the matter radius of 1.769 (1.524) fm for 3 He ( 3 H). 
In Ref. Jl2[ the effective cross section, cr e ff, was used to describe the interaction between 
the hadronic components in the virtual photon and the nuclear target. Here we take two 
models for cr e ff' the first (case 1) from Frankfurt and Strikman ||20|| , and the second (case 



2) the two-phase model of Ref. At large x(> 0.2) we used the structure functions of 



3 He and 3 H, obtained as a solution of the Faddeev equations for three body system P"T|.|2"2" 



Since the contribution of the nuclear binding and Fermi motion effects to the non-singlet 
combinations of the structure functions of 3 He and 3 H is small, we make an approximation 



that the binding and Fermi motion effect of Refs. [|Tl|,0l an d the pion cloud effect, discussed 
in the present work, are not correlated and, hence, contribute additively to 5F A=3 . 

We present our main results in Figs. and |] (the calculations were performed at Q 2 = 4 
GeV 2 ). To treat parton densities in the free proton and neutron realistically we used the 



CTEQ5L parametrization |23[. Since we know that the Coulomb effect is small in the region 
x > 10~ 4 , we set (3 = and show only the dependence of the flavor non-singlet structure 
function on a. 

As expected, the change of the pion fields leads to a considerable suppression of the non- 
singlet structure function in a nucleus. However, at very small x (typically ~ 10~ 3 — 10~ 4 ) 
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the reduction is not large compared with the enhancement caused by shadowing. (In the 
figures we present our results only for the range x G [1(T 2 , 0.8] in order to see the difference 
among the curves clearly. Note that the non-singlet structure function (divided by x) is of 
order 200 — 300 at x = 10~ 4 .) When the shadowing is turned off the difference of the proton 
and neutron structure functions gives the main contribution to the non-singlet structure 
function of the nucleus. In this case the change of the pion fields gives a sizable contribution 
at very small x (e.g., about 20% at x — 10~ 3 ). On the contrary, when the shadowing is 
switched on the main contribution to the non-singlet structure function at small x is given 
by a term proportional to (F% + F£) x (fs He — /s#), which is much larger than the effect of 
the change of the pion fields - the pion effect is at most 5% at x — 10 -3 . 

We can estimate the nuclear Gottfried integral, I(? e ' H (10~ 4 ), defined by Eq.fl3"3"|). If a 
is set to be (no change in the pion fields) we find I G He ' H (10~ 4 ) = 0.2953(0.3395) for case 
1 (2) - note that the CTEQ5L fit gives /^(HT 4 ) = 0.2403. When a = 0.1(0.2) we obtain 
4 He ' 3H (10" 4 ) = 0.2699(0.2444) for case 1, while it is 0.3142(0.2829) for case 2. Therefore, 
we expect that the change of the pion fields in the A = 3 system might lead to a reduction 
of the Gottfried integral by about 10%, compared with the value for the case where a = 0. 

We here give some comments: 

(1) Using the scaling assumption the pion distributions are shifted in the present calculation. 
The momentum fraction carried by pions is thus different from that in the free nucleon. For 
example, if < (3 < 1, we find (yf^ /3He (y)) - (1 + P)(vfoN(v))- Thus, in 3 He ( 3 H) the 
ratio of the momentum fraction carried by pions to that in the case where (3 = increases by 
about P/9 (4/?/9). The change of the pion momentum fraction would lead to changes of the 
momentum fractions carried by the nucleons and other mesons. In such a case one would 
have to construct a model where the momentum fractions are balanced. Note however, that 
when we set /3 — there is no change of the pion momentum fraction in a nucleus by virtue 
of the scaling assumption and pion number conservation. 

(2) We have treated the change of the pion fields in the nuclear medium semi-quantitatively. 
In principle, such a modification should be attributed to various correlation phenomena in 
a nucleus like Pauli blocking, exchange currents, short-range correlations etc. Melnitchouk 
and Thomas [^1 have reanalysed the nuclear shadowing effect on the deuteron structure 
function, including meson exchange currents, to extract the neutron structure function, 
and studied the Gottfried sum rule for the free nucleon. It should be possible to do such 
calculations for the three body system in the future. For a pair of mirror nuclei larger 
than the three body system, it would become important to consider Fock states including 
multi-pions, A isobars, other mesons and so on. That is also a very intriguing problem. 

In summary, we have estimated the effect of the medium modification of the pion fields 
on the flavor non-singlet structure function of the lightest mirror nuclei. We have found 
that the change of the pion fields produces a considerable suppression of the non-singlet 
structure function, and that the Gottfried integral is correspondingly reduced. In general, 
charge symmetry is broken and the Gottfried integral is divergent However, the in- 



dependence of the flavor non-singlet structure function of a pair of mirror nuclei would 
provide significant information on phenomena involving non-pQCD dynamics (such as the 
pion cloud) in the nuclear medium. Experiments on deep-inelastic scattering off various 



mirror nuclei should be possible in the future [25]. If we could vary the atomic number A 



and the difference between the proton and neutron numbers Y independently in measuring 
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the nuclear structure functions, it would stimulate a great deal of work which should lead 
to new information on the dynamics of nuclear systems. 

We would like to thank F. Bissey for providing us with the results for the structure 
functions of 3 He and 3 H in the large Bjorken x region. This work was supported by the 
Australian Research Council and Adelaide University. 
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FIG. 2. 5F^= 3 /x vs x (a = 0.1). The dotted, solid and dot-dashed curves show the results 
with = 0,0.1,0.2, respectively. 
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FIG. 3. Non-singlet structure function (divided by x) for case 1 discussed in the text. The 
dotted curve shows the non-singlet structure function for the free nucleon, while the dashed curve 
presents the result for the A = 3 system without any change of the pion fields. The upper (lower) 
solid curve is for the full calculation with a = 0.1(0.2). 
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